Medium— induced QCD cascade: democratic branching and wave turbulence 
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We study the average properties of the cascade of gluons that is generated by an energetic parton 
propagating through a quark-gluon plasma. We focus on the soft, medium-induced, emissions which 
control the energy transport at large angles with respect to the leading parton. We show that the 
effect of multiple branchings are important. In contrast to what happens in a usual QCD cascade 
in vacuum, medium-induced branchings are quasi-democratic, with offspring gluons carrying sizable 
fractions of the energy of their parent gluon. This results in a new mechanism for the transport of 
energy towards the medium, which is akin to wave turbulence with a scaling spectrum ~ 1/ ^/uj. We 
argue that the turbulent flow may be responsible for the excess energy carried by very soft quanta, 
as revealed by the analysis of the di-jet asymmetry observed in Pb-Pb collisions at the LHC. 
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One important phenomenon discovered recently in 
heavy ion experiments at the LHC is that of di-jet asym- 
metry^ a strong imbalance between the energies of two 
back-to-back jets. This asymmetry is commonly at- 
tributed to the effect of the interactions of one of the two 
jets with the hot QCD matter that it traverses, while 
the other leaves the system unaffected. Originally iden- 
tified [U [2] as missing energy, this phenomenon has been 
subsequently shown [3] to consist in the transport of a siz- 
able part of the jet energy by soft particles towards large 
angles. Some of the features of in-medium jet propaga- 
tion are well accounted for by the BDMPSZ mechanism 
for medium-induced radiation (from Baier, Dokshitzer, 
Mueller, Peigne, Schiff [4] and Zakharov (5). However, 
most studies within this approach have focused on the 
energy lost by the leading particle, while the LHC data 
call for a more thorough analysis of the jet shape for 
which the effects of multiple branching at large angles 
are important. Within that context, an important step 
was achieved in Ref. [6], where it was shown that, in a 
leading order approximation, one could consider succes- 
sive gluon emissions as independent of each other. This 
allows one to treat multiple emissions as a probabilis- 
tic branching process, in which the BDMPSZ spectrum 
plays the role of the elementary branching rate (THI] . 

Specifically, the differential probability per unit time 
and per unit z for a gluon with energy u to split into 
two gluons with energy fractions respectively z and 1 — z 
reads 
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where Pg^g{z) = [^hz) + ^(1 - z) - 2 
leading order gluon-gluon splitting function, (jcff = 
(J [1 — z(l — z)], with q the jet quenching parameter (the 
rate for transverse momentum broadening via interac- 
tions in the medium), and t^^_^{z,uj) is the time scale of 
the branching process. Note that we use light-cone (LC) 
coordinates and momenta, with the longitudinal axis de- 



fined by the direction of motion of the leading particle. 
Correspondingly, the 'energy' w truly refers to the LC 
longitudinal momentum p'^ and t to the LC 'time' x^. 
Eq. ([T]) apphes as long as ^ ^ t^^(z,w) < L, where L is 
the length of the medium, and (. is the mean free path be- 
tween successive collisions. The second inequality above 
implies an upper limit on the average energy of the off- 
spring gluons: zil — z)uj < Wc, where ojc = qL'^ /'^ is 
the maximum energy that can be taken away by a single 
gluon. It follows from Eq. ([T]) that the probability for 
having just one emission throughout the medium is (for 
z not too close to 1) ~ ay^uJc/ zui, where a = agNc/n. 
When this becomes of 0{1), i.e, when zuj < ojs = a^uJc, 
multiple branchings become important. 

It will be useful to express the energy of a radiated 
gluon in terms of its energy fraction x = uj/E and to 
replace the light-cone time t by the dimensionless variable 
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where Xc = iOc/E. We restrict ourselves here to the case 
E < LOc, i.e. Xc > 1, leaving the discussion of the E > ujc 
case to a forthcoming publication. Note that the maxi- 
mal value of r is r^ax — oi\f^c, corresponding to t — L. 
Then the branching probability ([T]) can be written as 
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where JC(z) = f{z)/[z{l - z)]^/^ = JC{1 - z) and f{z) = 

[l-z{l-z)f\ 

In this letter, we focus on one observable that charac- 
terizes the average properties of the in-medium cascade, 
namely the gluon spectrum, D{x,t) = x^. Within the 
probabilistic picture for multiple branchings, one easily 
derives the following evolution equation for D{x,t) 



dD{x,T) 



Az IC{z) 



D x,t] 



(4) 



2 



The initial condition corresponds to a single gluon car- 
rying all the energy, that is, D{x, r = 0) = 6{x — 1). We 
shall refer to the right hand side of Eq. Q as the "colli- 
sion term", and denote it as I[D]. Its physical interpre- 
tation is clear: The first contribution, which is non-local 
in X (except when x is close to 1), is a gain term: it 
describes the rise in the number of giuons at x due to 
emissions from giuons at larger x. Note that the func- 
tion D{x, t) has support only for < x < 1, which limits 
the first z-integral in Eq. ([4]) to a; < z < 1. The second 
contribution to the collision term, local in x, represents 
a loss term, describing the reduction in the number of 
giuons at X due to their decay into giuons with smaller 
X. Taken separately, the gain term and the loss term in 
Eq. Q have endpoint singularities at z = 1, but these 
singularities exactly cancel in between the two terms and 
the overall equation is well defined. 

For T ^ 1, we may attempt to solve Eq. Q in per- 
turbation theory, i.e., by iterations. Thus, by substi- 
tuting, in the collision term, D{x, r) by its initial value 
D^°\x) — 6{x — 1), one obtains (for a; < 1) 



D'^^\x,X) = 
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For T = Tmax, this is just the BDMPSZ spectrum. For 
reason that will become clear shortly, we refer to the 
small X part of this spectrum as the "scaling spectrum" , 
i.e., Dsc{x) ~ A priori, because one expects the 

small- a:; region of the spectrum to be populated by mul- 
tiple branchings, one would expect perturbation theory 
to break down when t > ^/x, and consequently the spec- 
trum to be strongly modified in this region. As we shall 
see, this is not at all the case: the scaling spectrum re- 
mains remarkably stable. In fact, if we were to continue 
solving Eq. Q by iteration, i.e., inserting the spectrum 
([5| into the collision term, one would find that after an 
obvious cancellation in the collision term, the remaining 
contribution is of the form of the scaling spectrum (see 
also [7] for a similar observation). 

The perturbative analysis turns out to be complicated 
by the non local character of the gain term (implying a 
non trivial coupling between the small- a; region and the 
region near a; = 1). In order to go beyond perturbation 
theory, and get insight into the non-perturbative features 
of Eq. (|4|, we have considered a simpler version of this 
equation, obtained by modifying the kernel to ]Co{z) = 
l/[z(l — z)]'^/^ (i.e, replacing the smooth function f{z) 
by 1 in Eq. ([s])). This simplification does not affect the 
singular behavior of the kernel near z = and z — 1, 
which determines the qualitative features of the solution, 
but it allows us to solve Eq. Q exactly, via a Laplace 
transform. The solution reads 
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The essential singularity at a; = 1 is a non perturbative 
effect that can be understood as a Sudakov suppression 



FIG. 1: Plot (in Log-Log scale) of •^Do{x,t), with D(i{x,T) 
given by Eq. ([gJi, as a function of x for various values of r 
(r = 0.01,0.02,0.1,0.2,0.4,0.6,0.9). [Color online] 



factor |5] (i.e. the vanishing of the probability to emit 
no gluon in an arbitrary small time). Aside from this 
exponential factor, one recognizes the scaling spectrum 
to which Dq{x,t) is proportional at small x. This is 
illustrated in Figjl] we see that the scaling spectrum 
is established early on (in agreement with our previous 
discussion), and remains stable as time progresses. For 
small times, its amplitude grows linearly with r: the sys- 
tem can then be viewed as a radiating source located at 
X < 1, feeding all the small-a; modes. As time passes 
the source weakens and eventually disappears into the 
left moving "shock wave' ' visible in Fig. [l] The time 
dependence of the solution is governed by the factor 
T exp[— 71x^/(1 — a;)] which has an a;-dependent maxi- 
mum: For X = 0.8 the maximum occurs at t ~ 0.18; 
this is the time where the peak starts to weaken ap- 
preciably, and the shock wave develops. For a: = 
the maximum occurs at t = l/^/2n ~ 0.4; after that 
time the spectrum decreases exponentially fast. The to- 
tal amount of energy stored into the modes is given by 
£o{t) = dxDQ{x,T) = e~^'^^ . It is nearly constant for 
^ l/v^i ^i^d decreases very rapidly after that. Be- 
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cause the spectrum keeps the same shape as one goes to 
smaller and smaller x, there is no accumulation of energy 
at any value a; > 0. Energy conservation then implies 
that a 'condensate' develops at a; = 0, where the excess 
energy coming from the large x region gets stored. We 
shall return later to the physical interpretation of this 
phenomenon. 

At this point, it is instructive to analyze an auxiliary 
problem — that of a system driven by a permanent source 
of energy localized at a; = 1. Consider then the equation 
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FIG. 2: The function ^/xDtb{x,T) (for A = 1) at early time, 
where it resembles Fig. [l] till late time when it approaches 
the steady state and saturates at the value A/2n at small x. 
The values of r are: 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1. 



solves this equation with initial condition Dti,{x,T = 
0) = 0. (Notice that the derivative of Dth{x,T) is equal 
to Dq(x,t), to within the multiplicative constant A.) 
By comparing Fig. [T] and Fig. [2] one sees that the be- 
haviors of D{x,t) and D^^{x,t) are remarkably similar 
at small r. However the most remarkable property of 
Dth{x, r) is that it converges to a time-independent func- 
tion, Dst{x) = {A/2n)/y/x{l — x). To understand this, 
we note that Dst{x) annihilates (exactly) the collision 
term, i.e., I[Dst] = 0, as can be verified by an explicit 
calculation. As time goes on, the solution Dti,{x,T) is 
gradually driven to Dst{x), but since this fixed point so- 
lution reduces to the scaling spectrum at small x, one 
observes very little change in this small- a; region, except 
for an overall time-dependent scaling. 

The stability of the scaling spectrum, confirmed by 
this simple calculation, hides important physics that is 
revealed by calculating the flow of energy that gets trans- 
mitted per unit time from the region x > xq to the region 
X < Xq. If we denote by £{xo, t) = j^^^ AxD{x, r) the to- 
tal energy that is contained in the modes with x > xq and 
recognize that the rate of change of £{xq, t) is due both 
to a possible source of strength A localized at x = 1, and 
to the flow V{xq,t) at xq, we get the general expression 
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An explicit calculation for D — Dth yields 
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where erfc(x) denotes the complementary error function. 
In order to analyze the physical content of this expres- 
sion, it is actually useful to rewrite the integral of the 
collision term in Eq. ^ in the following form 

7'(xo,r)- / dz/C(z) / dx ^ (11) 
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At small times, ttt^ ^ 1, and for xq not too close to 
either or 1, one can use the expansion erfc(x) ~ 1 — 
2x/y/^ in Eq. ([lo| and get ^(xo, t) ~ 2At ^/aVll^^a^, 
which vanishes at Xq = 0: At early times, the energy 
provided by the source is used to populate the vari- 
ous modes (according to the BDMPSZ spectrum), and 
does not flow across the entire system. Note that this 
contribution is obtained from Eq. ( |11[ ) by substituting 
D{x,t) by At6{1 — x), as appropriate at small times. 
This yields indeed 7'(xo,t) 

2At^xo/{1 
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xo). As this last calculation reveals, this 
early time contribution to the flow involves, when xq <C 1 
very asymmetric branching (z < xq). For this reason we 
interpret it as "direct radiation" . 

As time goes on however, the distribution of energy 
among the various modes is such that gain and loss 
terms equilibrate locally, at which point a steady state is 
reached with all the energy provided by the source flow- 
ing throughout the entire system and leaving the pop- 
ulation of the various modes constant. The constant 
(xo-independent) contribution to the flow is easily ob- 
tained from Eq. (10) by evaluating V{xq,t) at xg = 0, 
i.e., Vixo 



0,t) - A{1 



-TTT This result can be 
recovered from Eq. (Ill with xq <C 1 by approximating 
D{x,t) with the scaling part of the spectriun ([s]), i.e., 

A{1 - c-''''^)/(2tt^/x), and noting that 
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What this second calculation demonstrates, is that, in 
contrast to what happens for the direct radiation, here 
the typical branchings involve the whole range of z values 



(about half of the value of the integral (12) comes from 



the range 0.15 ^ z < 0.85). We refer to this property as 
"quasi-democratic branching" . 
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FIG. 3: The flow 'P{xo, r) as a function of xq from early times 
to the time where a stationary state is reached {A = 1). The 
values of T are: 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1. [Color online] 

The properties that we have discussed, the existence of 
a steady solution when the system is coupled to a source, 
the scaling spectrum whose origin can be traced back to 
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energy conservation and the particular z-dependence of 
the kernel IC{z), the presence of a component of the flow 
that is independent of the position, the quasi-democratic 
branching, all these are properties that are very reminis- 
cent of wave turbulence [TD]. Of course, in the present 
case, the locality of interactions in momentum space, of- 
ten a crucial ingredient for establishing the existence of 
steady state solutions, is only marginally satisfied (via 
what we called quasi-democratic branching) . And indeed 
the separation of the energy flow between radiation and 
quasi-democratic branching, that we introduced above, 
is not a clear-cut one; in general both components are 
mixed in a complicated way. Still these two components 
are clearly visible in Fig. |3j where we recognize the typ- 
ical distribution of modes fixed by early time radiation, 
the nearly uniform flow that develops at small times, and 
the approach to the steady state at later times. 

These properties can be contrasted with those of a typ- 
ical QCD cascade described by the DGLAP equation. 
Such an equation can be recovered from Eq. Q with 
the substitution K,(z)ly/x \l\z(\ — z)], where now 
r = Q;ln(Q^/QQ), with the virtuality. The effect of 
the evolution is again to redistribute the energy initially 
concentrated at a; = 1 towards smaller x. However, in 
contrast to the case studied in this paper, all the energy 
of the cascade remains in the spectrum, there is no tur- 
bulent flow. In fact it is not difficult to estimate that, 
in this case, the dominant contribution to the flow is 
V{xo) xq ln(l/a;o) which vanishes when xq — > 0. 

Many of the features that we have uncovered by study- 
ing the source problem remain valid without the source, 
and for the general kernel. We return now to this initial 
setting. Limiting ourselves to small times, one finds, by 



calculating the flow from Eq. (11 1, that it takes the form 
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with V — 4.96 obtained from Eq. (12) with the full f{z 



in the numerator. One recognizes the two components 
that we discussed earlier, that is, the direct radiation and 
the turbulent flow. Although parametrically subleading 
(since proportional to r), the turbulent flow dominates 
over the direct radiation when xq < r^, that is, in the 
region where multiple branchings are known to be im- 
portant. The total energy transported by this turbulent 
flow can be estimated by integrating the second term of 



Eq. ( 13 ) over time. One gets 
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This energy, which is not carried by the particles present 
in the spectrum, is what we identified earlier as the en- 
ergy stored in a 'condensate' at a: = 0. In more physical 
terms, we associate this energy with that transferred to 
the medium in the form of very soft quanta. 



It is beyond the scope of this letter to present a de- 
tailed comparison with the data. However the follow- 
ing order-of-magnitude estimates should confirm the rel- 
evance of the present discussion for the di-jet asymme- 
try observed at the LHC. Using the conservative esti- 
mate Wc — 60 GeV (corresponding to q — IGeV^/fm 
and L ~ 5 fm), together with 6? ~ 0.1, one finds 
£flow — 15 GeV, a value that compares well with the 
observations. Indeed, the detailed analysis by CMS [3] 
shows that the energy imbalance between the leading and 
the subleading jet is compensated by an excess of semi- 
hard {pt < 8 GeV) quanta propagating at large angles, 
outside the cone defining the subleading jet. For the most 
asymmetric events, the total energy in excess is about 
25 GeV. Remarkably, most of this energy (about 75%) is 
carried by very soft quanta with pr < 2 GeV. This obser- 
vation would be difficult to reconcile with the hypothesis 
that these particles come from gluons in a BDMPSZ spec- 
trum (which would imply that most of the excess energy 
would be carried by the hardest gluons with energies < 8 
GeV) . But it could be naturally explained by associating 
these soft particles with those transported by the turbu- 
lent flow that we have discussed in this paper. 
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